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SUMMARY 

The  objective  of  the  research  is  to  obtain  data  that  can  be  used  to 
predict  the  concentration  of  a  passive  tracer  released  into  the  wake  of  a 
submarine  travelling  through  a  stratified  body  of  water.  To  develop  this 
prediction  capability,  information  is  required  on  the  wake  shape  and  size. 
Experiments  have  been  performed  in  which  a  grid  is  oscillated  in  a  stably- 
stratified  flow  to  produce  the  steady-state  counterpart  of  the  momentumless 
wake  of  a  self-propelled  submarine  vehicle.  A  pH  sensitive  indicator  was 
used  to  produce  a  neutrally  buoyant  tracer  to  visualize  wake  development  and 
subsequent  vertical  collapse  as  the  wake  turbulence  is  overcome  by  the  buoy¬ 
ancy  gradient.  Measured  velocity  profiles  in  the  simulated  wake  indicate 
that  it  is  nearly  momentumless  and  vertical  temperature  surveys  reveal  the 
degree  of  mixing  in  the  wake. 

The  wake  growth  before  and  after  collapse  and  the  distance  to  collapse 
have  been  correlated  by  using  power  laws  previously  applied  for  the  wakes  of 
re-entry  vehicles  and  a  theoretical  analysis  of  marine  wake  collapse  developed 
under  the  present  program.  The  theory  is  based  on  the  balance,  at  wake 
collapse,  between  the  potential  energy  of  the  mixed  fluid  and  the  energy  of  the 
turbulent  velocity  fluctuations  in  the  wake.  Experimental  data  available  from 
the  literature,  as  well  as  that  obtained  in  the  present  program,  have  been  used 
in  the  correlation.  In  its  most  fundamental  form,  the  correlation  gives  the 
ratio  of  the  vertical  extent  of  the  wake  in  stratified  flow  to  that  in  unstratified 
flow  as  a  functior  of  the  ratio  of  the  time  in  the  wake  (measured  from  the 
body)  to  the  Brunt- Vaisala  period. 

The  scaling  relations  established  for  predicting  stratified  flow  wake 
dimensions  reveal  that  the  important  parameters  are  the  Froude  number, 
defined  as  the  product  of  the  submarine  velocity  and  the  Brunt- Vaisala  period 
of  the  ocean  divided  by  the  initial  wake  size,  and  the  ratio  of  the  residence 
time  in  the  wake  to  the  Brunt- Vaisala  period.  Through  the  use  of  these 
parameters,  two  unique  curves  are  obtained  for  estimating  the  horizontal 
width  and  vertical  height  of  a  wake  in  a  stratified  flow.  Although  enhance¬ 
ment  of  tv>e  horizontal  wake  growth  due  to  vertical  collapse  has  been  observed, 
the  area  of  the  wake  is  typically  l«.ss  than  one  half  that  for  unstratified  flow. 
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From  the  measurements  and  analysis,  equations  have  been  developed 
for  predicting  submarine  wake  dimensions  as  a  function  of  submarine  velocity, 
initial  wake  size,  ambient  stratification,  and  time  after  submarine  passage. 

The  estimated  submarine  wake  dimensions  are  of  the  correct  magnitude 
although  the  equations  should  brt  regarded  as  preliminary  since  further 
investigation  is  required  m  several  areas. 

The  very  flat  and  wide  wakes  observed  due  to  stratification  lead  to 
the  estimate  that  the  concentration  of  a  passive  tracer  introduced  into  a 
submarine  wake  would  remain  high  for  long  times  after  submarine  passage. 
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1. 


INTRODUCTION 


The  wake  of  a  submarine  travelling  through  a  stably- stratified  ocean 
exhibits  a  behavior  considerably  different  from  the  more  familiar  wake 
of  a  vehicle  moving  in  an  unstratified  environment.  In  both  cases,  the  wakes 
grow  initially  with  a  circular  cross-section.  However,  in  the  stratified  flow 
case,  mixing  in  the  near  wake  destroys  the  ambient  temperature  gradient  so 
that  further  downstream  buoyancy  forces  begin  to  inhibit  the  vertical  growth 
of  the  wake.  Sometime  after  generation,  the  wake,  which  has  a  density 
gradient  different  from  the  surroundings,  collapses  vertically  as  the  fluid 
seeks  its  own  density  level.  The  decrease  in  vertical  extent  of  the  wake  en¬ 
hances  the  horizontal  growth  rate  so  that  eventually  a  very  thin  and  wide  wake 
is  produced.  This  type  of  wake  history  in  stratified  flows  has  generated 
interest  from  the  point  of  view  of  studying  possible  internal  waves  generated 
by  wake  collapse  as  well  as  tracking  passive  trace  elements  released  by  the 
submarine . 

The  objective  of  the  present  research  program  is  to  obtain  analytical 
and  experimental  data  that  can  be  used  for  predicting  the  concentration  (as  a 
function  of  downstream  distance)  of  a  passive  tracer  released  into  the  wake 
of  a  self-propelled  body  travelling  through  a  stably-stratified  medium.  To 
achieve  this  prediction  capability,  information  is  required  on  the  wake  shape 
and  size.  The  experimental  approach  utilizes  a  stratified  flow  in  which  a 
composite  grid  is  oscillated  to  produce  a  steady-state  counterpart  of  the 
momentumless  wake  of  a  self-propelled  body.  A  pH  sensitive  indicator,  in 
which  local  color  changes  are  generated  at  the  grid  by  electrical  impulses,  is 
used  for  direct  observation  of  the  wake  development,  diffusion  and  subsequent 
vertical  collapse  as  the  stable  stratification  counteracts  the  wake  turbulence. 

From  the  present  measurements  and  analytical  efforts  along  with  those 
of  other  investigators,  techniques  are  established  for  correlating  the  experi¬ 
mental  data  available  on  wake  collapse  in  laboratory  environments.  The 
scaling  laws  produced  in  this  manner  th<  a  are  used  to  give  preliminary 
estimates  of  the  growth  and  collapse  of  the  wake  behind  a  submarine  travelling 
through  a  stratified  ocean  for  comparison  with  full-scale  data. 


1 


2. 


EXPERIMENTAL  TECHNIQUES 


A  schematic  arrangement  of  the  flow  system  used  for  the  experiments 
is  given  in  F  gure  1.  The  tests  were  performed  in  a  water  tank  (210  x  10  x 
10  cm)  in  \  .ich  a  stratified  flow  with  a  very  nearly  linear  temperature  gradient 
is  generated  by  the  mixing  of  streams  of  water  at  different  temperatures.  The 
degree  of  stratification  is  determined  by  continually  measuring  the  vertical 
temperature  profile  through  the  use  of  small  thermistor  beads  that  traverse 
the  depth  of  the  water  at  selected  locations  at  regular  intervals. 

The  velocity  of  the  flow  is  measured  by  tracers  that  are  neutrally 
buoyant  at  every  point  in  the  fluid.  This  is  necessary  because  the  motion  of 
a  constant-density  dye  would  be  influenced  by  the  density  gradient,  giving  a 
distorted  velocity  profile.  The  technique  *  is  based  on  the  color  change  of 
an  indicator  solution  when  the  hydrogen  ion  concentration  (pH)  is  changed  at 
a  test  point.  The  solution  used  consists  of  water  containing  0.01%  thymol 
blue.  It  is  titrated  to  the  end  point  by  adding  a  small  amount  of  hydrochloric 
acid  and  is  initially  bright  orange  in  color.  Fine  tungsten  wires  stretched 
vertically  across  the  water  act  as  cathodes,  and  the  anode  is  formed  by  a 
stainless  steel  wire  along  the  tank  bottom.  Application  of  brief  voltage  pulses 
(90  volts)  to  the  two  electrodes  reduces  the  hydrogen-ion  concentration  at  the 
surface  of  the  tungsten  ,vire  corresponding  to  a  local  increase  in  pH  .  Columns 
of  dark  tracer  fluid  are  formed  at  the  surface  of  the  tungsten  wire  and  move 
with  the  fluid  to  yield  vivid  velocity  profiles.  Typical  flow  speeds  are  about 
one  cm/ sec. 

As  shown  in  Figure  1,  a  wake  analogous  to  that  of  a  self-propelled  body 
is  generated  in  the  water  channel  test  section  through  the  use  of  an  oscillating 
grid  approximately  one  quarter  inch  in  diameter.  The  grid  is  pulsed  electric¬ 
ally  to  produce  a  local  pH  -  color  change  which  traces  the  wake  growth  down¬ 
stream  of  the  grid.  Measurements  of  the  wake  dimensions  in  both  the  horizontal 
and  vertical  planes  are  made  at  different  distances  from  the  grid  for  both  un¬ 
stratified  and  stratified  flow.  In  the  unstratified  case,  the  wake  cross-section 
at  any  point  is  circular.  Whe.'  flow  is  stratified,  buoyancy  causes  the  wake, 
which  initially  grows  at  the  same  rate  in  all  directions,  to  collapse  in  the 
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vertical  direction  and  flatten  out  horizontally  producing  a  very  wide  and  thin 
elliptical  or  rectangular  shape. 

The  determination  of  the  wake  dimensions  as  a  function  of  wake  position 
and  ambient  stratification  is  vital  for  the  prediction  of  the  concentration  of  a 
passive  tracer  released  at  the  body.  The  research  discussed  in  this  rcsport 
has  been  directed  towards  this  objective. 

3.  WAKE  GROWTH  IN  UNSTRATIFIED  FLOW 

The  growth  of  a  wake  downstream  of  a  body  travelling  through  a 
medium  in  which  there  is  no  temperature  gradient  has  been  well  documented 
in  the  literature.  Theoretical  analyses  along  with  low-speed  tests  in  water 
and  hypersonic  tests  in  ballistic  ranges  and  wind  tunnels  have  established 
relations  for  wake  growth  of  the  form 


where  'S'm  is  the  wake  width  (circular),  DB  the  diameter  of  the  body,  CD 
the  drag  coefficient  of  the  body  and  X  the  distance  in  the  wake  from  the  body. 
For  unpropelled  bodies,  measurements  of  the  turbulent  wake  growth  rates 
yield  values  of  about  unity  for  K  and  1/3  for  the  exponent  n  ,  in  good  agree¬ 
ment  with  theoretical  predictions . 

For  the  turbulent  wakes  of  bodies  with  hydrodynamical  self -propuls ion, 
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theoretical  analytes  yield  values  for  n  of  1/4  for  plane  flow  and  1/5  for 
axially  symmetric  v aw.  At  large  distances  behind  the  body,  the  wake  may 

4 

asymptote  to  a  constant  size  ,  although  this  has  not  been  verified  experimentally. 

5  .9 

Various  attempts  to  model  a  submarine  wake  have  been  reported  in 
which  a  two-dimensional  unsteady  flow  is  produced  through  the  use  of  paddle- 
type  mixers  in  a  tank  containing  quiescent  water.  In  these  instances,  the 
time  t  after  mixing  is  equivalent  to  the  distance  in  the  wake  X  in  Equation  (1). 

It  has  been  observed  on  the  present  research  program  that  these  measurements 
can  be  plotted  in  the  form  /  D  vs.  t/  D  for  the  cases  in  which  the  water 
is  unstratified  to  give  good  correlation  of  the  data,  with  being  propor¬ 

tional  to  (t/  D )  /3  ,  where  D  is  the  initial  diameter  of  the  mixed  region. 
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The  value  of  1/3  for  the  exponent  is  highland  indicates  that  the  mixed  regions  1 
generated  in  the  two-dimensional  unsteady  experiments  are  not  strictly 

l 

equivalent  to  those  at  large  distances  behind  a  self-propelled  body,  the  wake 
of  which  would  gr’ow  at  a  much  slower  rate. 

,  In  Figure  2,  it  is  shown  thalt  tests  reported  by  Schooley  and  Stewart^ 

*«  I 

using  a  self-propelled  model  yield  a  wake  in  unstratified  water  that  can  be 

correlated,  by  a  relation  =  l.O(-g-)  obtained  for  hypersonic  spheres.  The 

high  value  of  1/3  for  the  exponent  indicates  that  the  wake  growth  in  this  cia.se 

is  determined  largely  by  the  effects  of  body  drag  and  the  propeller.  Similarly, 

4  ■ 

experiments  by  Naudascher  ,  utilizing  a  jet-disc  model  in  air  to  produce  a 
steady-state  counterpart  of  the  momentumless  wake  of  a  self-propelled  body, 
give  an  initial  growth  rate  of  ■—  =  j.5(-g-)  .  This  very  rapid  growth  may  be 

due  to  the  influence  of  the  jet  for  which  the  exponent  n  has  a  value  of  unity. 

i 

For  distances  in  the  wake  greater  than  about  20  diameters,  the  wake  growth  in 
Naudascher's  experiment  can  be  represented  by:  &  7.2  (—)  consistent  with 

the  growth  rate  anticipated  for  a  momentumless  wake. 


In  addition  to  the  interpretation  of  previous  work  on  marine  wake  Simula 
tiou  in  terms  of  porter  law  growth1  rates,  measurements  have  been  made  of 

i  , 

the  wake  growth  downstream  of  an  oscillating  grid  in  a  one  cm/ sec  flow  of  un¬ 
stratified  water.  /The  wake  growth  observed  using  the  pH-color  change  tracer 
is  plotted  in  Figure  2  frpm  which  it  is  observed  that  the  >vake  diameter’  can  be 

i*.  '  -  *  * 
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expressed  by  the  i  elation  w  I  3  (  ~)  ' .  The  wake  diameter  reaches  a  size 

of  roughly  four  times  the  initial  diameter  at  one  hundred  diameters  downstream. 
The  measurements  indicate  that  the  wake  may  reach  an  asymptotic  size  al¬ 
though  furtker  tests  are  required  to  clarify  this  point.  The  vc.  <e  1/4  for  the 
exponent  in  the  growth  law  is  closer  to  the  value  anticipated  for  the  *’ake 
growth  rate  behind  a  self-propelled  vehicle  than  the  value  1/3  observed  i. 
other  simulation  experiments.  Waj^e  velocity  profile  measurements  obtained 
using  the  neutrally-buoyant  thymol  blue  tracer  indicate  little  if  any  velocity 
defect  confirming  that  the  grid  wake  is  nearly  momentumless  on  the  average. 


Computer  results  from  the  theoretical  analysis  of  Ko**  are  reported 
in  the  form  of  several  different  curves  for  Froude  numbers,  Fr  of  0.5,  1.0 
and  10.0.  It  is  shown  in  Figure  2  that  these  can  be  expressed  by  one  curve 
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of  the  form  =  1- 25  (-g- )°'IZ  during  the  early  growt1  period  before  the 

effects  of  stratification  are  felt.  This  relation  compares  very  well  with  the 

X  O  2.5  '  i 

equation  -=*■  »  f.3  (-g-)  observed  in  the  present  experiments  for  unstratified 
flow . 


The  energy  of  the  turbulence  generated  in  the  wake  by  the  oscillating 

grid  can  be  estimated  by  using  the  •  -verified^  result  cLlr„/dt  ,  where 

'  .  1  '  ' 

U1  is  a  characteristic  turbu1  .t  velocity  fluctuation,  and  is  the  width 

fl 

of  the  wake.  For  the  growth  rate  observed  in  the  experiments,  i.e.  —■  - 

'  v  ■  54  D  " 

this  yields  -  0-325  (•—)  for  the  rate  at  which  the  turbulence  intensity 

decays  downstream  of  the  oscillating  grid.  ,  In  Figure  3  this  relation  for  u 

4 

is  shown  to  agree1  remarkably  well  with  Naudascher's  hot  wire  anemometer 
measurements  in  air  of  the  turbulence  intensities  u  ,  v'  ,  uj'  on  the  axis  behind 
a  jet-disc  model  which  is  used  to  simulate  the  momentumless  wake  of  a  self- 
propelled  body.  This  agreement  confirms  the  hypothesis  that  the  energy  of  th^ 

■  1  ,  I 

turbulence  at  different  positions  in  the  wake  can  be  estimated  from  measure¬ 
ments  of  the  rate  at  which  the  wake  dimensions  change  with  time,  and  that  the 
turbulence  decay  rate  can  be  expressed  by  a  straightforward  power  law. 


The  results  obtained  from  the  laboratory  measurements  of  wake  growth 
and  turbulence  intensity^  shown  in  Figures  2  and  3,  can  be  used  to  estimate 

I 

the  size  and  energy  of  the  wake  of  a  submarine  travelling  through  an  unstratified 

|  I 

ocean.  The  predicted  diameter  of  the  wake  behind  a  typical  vehicle  would  in¬ 
crease  from  150  feet  at  two  miles  doWnstream  to  225  feet  at  ten  miles.  In  the 
same  distance,  the  turbulent  velocity  fluctuations  in  the  wake  of  a  vehicle 
travelling  at  10  knots  would  decay  from  approximately  2.5  cm/ sec  at  two  miles 
to  0.5  cm/sec  at  10  miles.  The  wake  energy  represented  by  the  square  of  the 
turbulence  velocity  plays  an  important  role  in  determining  the  collapse  of  the 
wakfe  when  the  flow  is  stratified., 


4.  ,  WAKE  GROWTH  AND  COLLAPSE  IN  STRATIFIED  FLOW 

4.1  Wake  Collapse  Time 

The  wake  of  a  body  trav'-  ,.ng  through  an  environment  that  is  stably 
stratified  can  be  profoundly  ir  lu^nced  by  the  density  gradient.  Immediately 
downstream  of  the  vehicle,  the  wake,  which  may  contain  fluid  of  nearly 
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constant  density,  will  grow  at  the  same  rate  in  all  directions,,  as  in  the  un¬ 
stratified  flow  case,  (Figure  2).  However,  as  the  turbulent  energy  of  the 
wake  decays  with  increasing  distance  from  the  body  (Figure  3),  the  restoring 
action  of  buoyancy  begins  to  inhibit  the  vertical  expansion  of  the  wake  and  at 
the  same  time  enhances  the  horizontal  growth.  At  some  point  behind  the 
body,  the  wake  reaches  a  maximum  vertical  size  followed  by  a  collapse  as 
the  fluid  returns  under  the  action  of  gravity  to  the  level  at  which  its  density 
is  the  same  as  the  environment. 
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Many  investigations  ,  primarily  two-dimensional  unsteady  ones , 
have  been  directed  toward  determining  the  distance  or  time  from  the  body  fo 
wake  collapse.  A  parameter  used  to  characterize  the  phenomenon  is  the 
Brunt-Vaisala  period  T  defined  by 
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where  3^  /  dz  is  the  ambient  density  gradient. 

Since  T  is  the  period  at  which  a  parcel  of  displaced  fluid  oscillates 

about  its  equilibrium -density  position,  it  would  be  reasonable  to  expect  that 

the  time  to  wake  collapse  tc  would  be  a  function  of  T  .  In  a  recent  review 
13 

by  Sundaram  ,  apparently  conflicting  results  were  presented  in  an  attempt 

to  relate  tc  to  T  .  However,  it  was  observed  on  the  present  research 

program  that  this  disagreement  between  results  from  different  investigators 

5  8 

can  be  attributed  to  the  fact  that  Van  de  Wat'  ring  ’  expresses  the  Brunt- 

6  7 

Vaisala  period  per  radian  and  Schooley  *  per  cycle.  In  addition  Schooley 
measures  the  time  to  collapse  tc  from  the  start  of  mixing  and  Van  de  Watering 
from  the  end  of  mixing.  As  shown  in  Figure  4a,  this  can  seriously  influence 
the  results  since  typical  values  for  T  are  between  5  and  24  sec/cycle  with  a 
typical  mixing  time  of  three  seconds. 

When  allowance  is  made  for  the  different  definitions  of  T  and  zero  time 
in  Figure  4  a  tc  /  T  from  the  two-dimensional  unsteady  experiments  is 
observed  to  lie  between  1/4  and  1/2  at  low  values  for  T  where  the  mixing 
time  is  equal  to  or  greater  than  tc  .  The  mixing  time  becomes  a  decreasing 
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fraction  of  the  collapse  time  with  increasing  T  ,  and  the  ratio  tc  /  T 
approaches  a  value  of  1/3  at  T  ■»  24  sec/cycle,  the  upper  limit  of  the  meas¬ 
urements.  For  the  experiments^  with  the  self-propelled  model  in  which  the 
mixing  time  is  very  small  relative  to  the  Brunt-Vaisala  period  T  (which 

itself  is  only  three  seconds),  ta  IT  in  Figure  4  a  equals  1/3.  Measurements 
14  15 

reported  ’  for  towed  plates  are  also  shown  to  give  approximately  the  same 
value  for  the  ratio  tc  /  T  in  Figure  4a. 

In  the  present  series  of  experiments,  the  wake  grows  downstream  of 
an  oscillating  grid  that  is  operated  continuously  in  a  flow  of  thermally- 
stratified  water  as  depicted  in  Figure  i.  Possible  influences  of  mixing  time 
are  eliminated  and  the  collapse  of  the  wake  is  not  constrained  by  the  presence 
of  the  mixer  in  the  center  as  in  the  two-dimensional  unsteady  experiments. 
Measurements  were  performed  for  Brunt-Vaisala  periods  between  8  and  60 
extending  the  range  in  T  closer  to  the  high  values  (in  the  hundreds)  measured 
in  the  ocean.  The  results  shown  in  Figure  4  b  indicate  that  the  collapse  time 
is  roughly  one-third  the  Brunt-Vaisala  period  over  the  complete  range  of  T  . 

It  should  be  emphasized  that  the  collapse  time  in  all  the  experiments  can  be 
estimated  only  approximately  since  the  arresting  of  the  vertical  wake  growth 
followed  by  the  collapse  is  a  gradual  process  spread  out  over  a  considerable 
period  of  time. 

4.2  Analysis  of  Wake  Collapse 

To  predict  the  concentration  of  a  passive  tracer  introduced  in  the  wake 
of  a  self-propelled  body  travelling  through  a  stratified  medium,  it  is  first 
necessary  to  establish  methods  for  predicting  the  size  of  the  wake.  As  part 
of  the  present  program  of  research,  a  theoretical  analysis  of  wake  collapse 
was  carried  out.  In  the  vicinity  of  collapse,  the  turbulent  energy  in  the  wake 
(which  makes  it  grow)  is  assumed  equal  to  the  potential  energy  (which  makes 
it  collapse) . 

For  a  wake  of  rectangular  cross-section  with  a  linear  density  gradient, 
the  potential  energy  can  be  approximated  by  hi  and  the  kinetic 

energy  by  -£v  /l  where  is  the  vertical  height  of  the  wake, 

is  the  horizontal  width,  9  is  the  ambient  density,  u.'  is  the 
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turbulent  velocity  fluctuation,  q  is  the  acceleration  due  to  gravity,  a=  -g- 

«  1  3o  i  Y  iZ  *' 

is  the  ambient  density  gradient  and  -y^-j  is  the  density  gradient  in  the 

wake  caused  by  the  mixing.  If  the  wake  were  completely  mixed,  the  wake 

would  have  a  constant  density  and  /3  would  equal  zero.  The  assumption  of  a 

rectangular  wake  will  be  shown  later  to  be  reasonable. 


The  representation  of  the  potential  and  kinetic  energies  by  the  above 
expressions  should  be  regarded  as  giving  a  pararnet  ic  rather  than  an  exact 
dependence.  Equating  the  two  approximate  equations  for  the  kinetic  and 
potential  energies  yields  the  following  relation  for  the  vertical  height  of  the  wake. 


Cu ') 

3  (*-/}) 


(3) 


From  measurements  {Figure  3),  it  has  been  established  that  for 
momentumless  wakes  in  unstratified  flow  the  turbulent  intensity  u/  can  be 
equated  to  the  wake  growth  rate  by  u!  «  •  Combining  this  relation 

with  Equation  (1),  and  taking  %  -  u-t  where  u.  is  the  vehicle  velocity, 
yields 


U. 


a 
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Combining  Equations  (4),  (3)  and  (2)  gives 


K 


0.  4- 0  rt 


(5) 


In  other  words,  the  ratio  of  the  vertical  width  "£v  of  the  wake  under 

stratified  flow  conditions  to  the  width  ^  when  the  flow  is  not  stratified  can 

be  expressed  as  a  simple  function  of  t  /T  (the  ratio  of  time  t  =  X./u 

between  the  body  and  the  wake  location  to  T  the  Brunt-Vaisala  period). 

$  */2 

The  degree  of  mixing  in  the  wake  enters  in  the  term  (1  -  -jj-} 

A  physical  understanding  of  the  t/T  scaling  can  be  obtained  by  noting 
tloc  interactions  involving  flow  stratification  are  often  characterized  by  the 
local  Richardson  number  R+  which  is  defined  in  Equation  (6). 
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*  =  (u)‘  e  3Z  (6) 

The  parameter  t  /  T  is  directly  related  to  RjL  since  Equation  (6)  leads  to 


R* 


(?) 


Essentially  the  Richardson  number,  R^  ,  represents  the  ratio  of 

the  potential  energy  to  the  turbulent  energy.  Using  Equation  (7)  in 

v 

Equation  (5)  reveals  that  $rvl4rn~  0.8U  /(Ri)  *0  -  fi/a.)  .  Hence,  as  physically 
anticipated,  an  increase  in  potential  energy  (higher  Richardson  number)  de¬ 
creases  the  vertical  extent  of  the  wake  while  an  increase  in  turbulence  (lower 
Richardson  number)  increases  the  vertical  wake  size. 


4.3  Measurements  of  Degree  of  Mixing  in  Wake 


An  important  parameter  revealed  in  the  analysis  of  watte  collapse. 
Section  4.2,  is  the  degree  of  mixing  in  the  wake,  ;/6  / OL  .  The  effect  of 
mixing  on  the  time  to  wake  collapse  ic  can  be  determined  by  taking 
in  Equation  (5).  This  yields 


te  o.4o  a 

T  ‘  -  &r  c) 

In  Figure  5,  the  ratio  tc  /  T  is  plotted  against  the  degree  of  mixing 
/3 /(X  for  A  equal  to  1/3  and  1/5.  For  the  completely  mixed  case, 

(  ft>/cc  =  0  ,  and  a  -  7/4  ),  tc  /  T  *  0.10  whereas  for  the  slightly 
mixed  case,  (/S/(X  =  0,1  ),  tc/T  0.32  .  This  is  the  value  observed 

in  Figure  4  that  best  correlates  the  available  experimental  data  on  wake 
collapse  time.  As  shown  in  Figure  5,  an  almost  identical  variation  in  tc/T 
with  /3  /  OC  has  been  observed  in  the  computer  results  from  Ko's^  theoreti¬ 
cal  analysis  in  which  the  value  of  ~  0.1  also  was  observed  to  best 

fit  the  measurements. 
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To  determine  the  degree  of  mixing  in  the  wake  as  a  function  of  down¬ 
stream  distance,  measurements  of  the  vertical  temperature  profiles  both 
inside  and  outside  of  the  wake  have  been  obtained  using  small  thermistor  beads 
that  traverse  the  flow.  The  results  are  presented  in  Figure  6  for  a  Brunt  - 
Vaisala  period  of  16.  The  ambient  profile  given  in  Figure  6  a  shows  a  temper- 
mature  gradient  of  roughly  3/4  0  C/cm.  In  Figure  6b,  a  temperature  probe  at 
a  fixed  vertical  position  reveals  temperature  fluctuations  in  the  wake  of  about 
0. 1°  C  compared  with  the  smooth  profile  obtained  in  the  undisturbed  stream 
when  the  grid  is  not  oscillated.  This  fluctuation  is  equivalent  to  roughly  10 
to  20%  of  the  ambient  temperature  difference  across  the  wake  at  this  point. 

The  probe  is  located  just  inside  the  upper  edge  of  the  wake  and  the  drop  in 
temperature  shown  in  Figure  6  b  when  the  grid  is  oscillated  confirms  that  the 
water  is  being  mixed.  This  can  be  seen  in  the  wake  temperature  profile  in 
Figure  6  c  in  which  a  wake  of  nearly  constant  temperature  has  been  generated 
near  the  grid.  As  shown  in  Figures  6d,  e,  and  f,  the  wake  temperature  profile 
returns  closer  to  the  ambient  one  with  increasing  distance  downstream  of  the 
grid  although  large  scale  motions  hinting  at  interval  waves  are  observed. 

From  measured  profiles  such  as  those  given  in  Figure  6,  the  degree 
of  mixing  fl/oi  in  the  wake  was  calculated.  The  results  are  presented  in 
Figure  7  as  a  function  of  time  after  wake  generation.  Although  considerable 
uncertainty  exists,  it  is  apparent  that  by  the  start  of  collapse  at  about  t/T  -  U3 
(Figure  4),  fl/Ct  is  approaching  a  value  not  much  less  than  unity.  This  is 
consistent  with  analytical  observations  discussed  in  Section  4.2  requiring  a 
value  of  fi/oL  =  0.9  to  fit  theory  to  experiment. 

4.4  Measurements  of  Wake  Collapse 

In  addition  to  the  analysis  of  marine  wake  collapse  and  the  correlation 

of  earlier  results  carried  out  on  the  present  contract,  measurements  of  wake 

growth  and  collapse  were  obtained  for  Brunt-Vaisala  periods  up  to  60.  This 

compares  with  a  value  of  3  for  the  self-propelled  model  experiments*^  and 
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6  to  24  for  the  two-dimensional  unsteady  experiments  that  are  influenced 
by  the  mixing  time  and  the  presence  of  the  mixer  at  the  center  of  the  turbulent 
region.  The  present  experiment  was  run  continuously  with  the  temperature 
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gradient  allowed  to  decrease  gradually  with  increasing  time  so  that  T 
varied  from  about  8  sec/ cycle  at  the  start  of  the  experiment  to  60  at  the  end, 
several  hours  later.  During  this  period,  the  wake  collapse  was  observed  to 
move  downstream  as  predicted  by  the  relation  tc  /  T  ^  7/3 

Measurements  of  the  vertical  and  horizontal  growth  of  the  wake,  ob¬ 
tained  using  the  pH-color  change  technique  and  the  experimental  configuration 
shown  in  Figure  1,  are  presented  in  Figure  8  for  a  Brunt-Vaisala  period  of  55. 
The  measurements  are  normalized  by  D  ,  t;  e  initial  wake  diameter  which  is 
approximately  equal  to  the  grid  diameter.  From  the  grid  to  an  %/D  of  15 
which  corresponds  to  a  value  for  t /T  of  0.15,  the  vertical  and  horizontal 
^  size  of  the  wake  are  approximately  equal  to  the  diameter  meas¬ 

ured  in  the  unstratified  flow  case,  Figure  2.  Further  downstream,  the  vertical 
size  of  the  wake  reaches  a  maximum  of  about  three  diameters  and  then  de¬ 
creases  as  the  wake  collapses.  This  causes  the  horizontal  extent  of  the 
wake  to  grow  at  a  faster  rate  than  observed  in  the  unstratified  flow  case. 

The  vertical  size  of  the  wake  appears  to  reach  an  asymptote  of  roughly 
1.7  times  the  initial  diameter  by  a  time  ce  about  one  Brunt-Vaisala  period. 

The  cross-sectional  area  of  the  wake  calculated  from  the  measurements 

in  Figure  8  is  shown  in  Figure  9  as  a  function  of  time  after  wake  generation. 

The  wake  area  is  normalized  by  the  initial  area  of  the  wake  .  For 

unstratified  flow,  the  circular  wake  grows  according  to  the  relation 

=  (<■  3)  (  )  /j  obtained  from  the  measurements  in  Figure  2.  For 

stratification,  an  assumption  must  be  made  concerning  the  shape  of  the  wake. 

Ko^  assumes  an  elliptical  wake.  In  Figure  9,  this  assumption  is  shown  to 

lead  to  an  unrealistic  decrease  in  wake  area  just  after  collapse  at  t/T  c  7/3  . 

The  same  behavior  is  observed  when  the  v/ake  area  is  calculated  for  meas- 
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urements  reported  by  other  investigators  .  A  rectangular  wake  will 
give  values  for  A  w  which  are  too  high  at  collapse  but  which  will  likely 
be  realistic  at  later  time. 

The  most  probable  variation  of  wake  area  with  time  is  shown  by  the 
dotted  line  in  Figure  9.  The  initially  circular  wake  starts  to  collapse  in  the 
form  of  an  ellipse  but  quickly  approaches  a  rectangular  cross-section. 


initial  growth  r?/.te  only  slightly  so  that  «*  0.9  to  1.0  before  collapse 

at  t/T  ^  1/3.  The  density  gradient  inside  the  wake  at  this  time  and  not 
that  at  t  =  0  determines  the  wake  collapse.  From  Figure  9»  it  is  seen  that  the 
wake  area  at  t/T  =  1/3  is  roughly  ten  times  the  initial  area.  At  i  =  0,  the 
wake  may  well  be  completely  mixed  with  y6  / oc  3  0  .  However,  with  growth, 
the' entrainment  of  such  a  large  quantity  of  ambient  fluid  into  the  wake  would 
be  expected  to  produce  a  density  gradient  at  collapse  not  much  different  from 
the  ambient  gradient,  as  observed.' 


To  estimate  the  vertical  extent  of  a  wake, during  the  collapse  phase, 
Equation  (5)  can  be  combined  with  the  measured  growth  rate  in  unstratified  flow 
P  *  K  (X/B)  to  yield 


ix  (A)  = 
D  '  F 


o  .40  Kn 
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1  D  ' 
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where 


u-T 


F  -  -jj-  is  a  Froude  number  defined  using  the  submarine  velocity 

,  and  the  initial  wake  diameter 


u,  ,  the  ambient  Brunt-Vaisala  period  I 
B  .  Comparison  with  Ko's^  definition  of  a  Froude  number  Fr  -  a.rT /2TCB 
where  U-y  ,  the  initial  turbulent  intensity  is  taken  by  Ko  to  be  equal  to 
0.25  ul  ,  gives  a  relation  between  the  two  Froude  numbers  of  F  «  25  Fr  . 

'  S’  i 

Following  Equation  (9),  (-=-)  is  plotted  against  %/B  in  Figure  11. 

11  Ur' 

For  comparison,  Ko's  theoretical  curves  from  wake  collapse  on  downstream 
for  Fr  of  0.5,  1.0  and  10  are  also  given  in  the  figure.  When  plotted  in  this 
fashion,  Ko's  results  for  different  Froude  numbers  fall  into  one  curve  given 
by  ~trv/B(F)  =  0.5 7(-q~)  *  For  the  present  analysis,  as  shown  in  Figure  11, 
this  corresponds  to  Equation  (9)  with  K  =  1.3,  /$/&  =  0.94,  and  n  =  1/5. 


A  simple  scaling  law  for  the  post  collapse  behavior  of  the  wake  is  thus 
obtained  from  the  two  analyses  when  the  vertical  size  of  the  wake  is  divided 
by  the  Froude  number,  F  ,  which  is  proportional  to  the  ratio  of  the  turbulence 
energy  to  the  potential  energy.  In  physical  terms,  at  a  given  wake  position 
X/D,  the  vertical  extent  of  the  wake,  ~$rv  / D  ,  which  is  directly  proportional 
to  the  Froude  number,  F  ,  increases  as  the  wake  turbulence  increases  with 
submarine  velocity  u,  and  decreases  as  the  potential  energy  (D/T)  increases. 
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initial  growth  r;>.te  only  slightly  so  that  /  ^<»  *»  0.9  to  1.0  before  collapse 

at  t/T  ^  1/3.  The  density  gradient  inside  the  wake  at  this  time  and  not 
that  at  t  =  0  determines  the  wake  collapse.  From  Figure  9*  it  is  seen  that  the 
wake  area  at  t/T  =  1/3  is  roughly  ten  times  the  initial  area.  At  t  =  0,  the 
wake  may  weU  be  completely  mixed  with  /S  / oc  *  0  .  However,  with  growth, 
the' entrainment  of  such  a  large  quantity  of  ambient  fluid  into  the  wake  would 
be  expect  jd  to  produce  a  density  gradient  at  collapse  not  much  different  from 
the  ambient  gradient,  as  observed.' 


To  estimate  the  vertical  extent  of  a  wake, during  the  collapse  phase, 
Equation  (5)  can  be  combined  with  the  measured  growth  rate  in  unstratified  flow 
«  K  CX/B)n  to  yield 


i*  (±)  = 

D  ^  F ' 


5.40  Kn 


n- 1 


o-4) 


* 


(9) 


where  F  = 


~~  is  a  Froude  number  defined  using  the  submarine  velocity 

u.  ,  the  ambient  Brunt-Vaisala  period  T  ,  and  the  initial  wake  diameter 

B  .  Comparison  with  Ko's ^  definition  of  a  Froude  number  Fr  =  U-rT /2TCJD 

»  ! 

where  (JL:r  ,  the  initial  turbulent  intensity  is  taken  by  Ko  to  be  equal  to 
0.25  ul  ,  gives  a  relation  between  the  two  Froude  numbers  of  F  »  25  Fr  . 

£  i 

Following  Equation  (9),  -pp  /-r-)  is  plotted  against  % /i?  in  Figure  11. 

11  Ur' 

For  comparison,  Ko's  theoretical  curves  from  wake  collapse  on  dow-stream 
for  fr  of  0.5,  1.0  and  10  are  also  given  in  the  figure.  When  plotted  in  this 
fashion,  Ko's  results  for  different  Froude  numbers  fall  into  one  curve  given 
by  £v/B(F)  =  0. 57(-£-)  .  For  the  present  analysis,  as  shown  in  Figure  11, 
this  corresponds  to  Equation  (9)  with  K  =  1.3,  /3/oi  =  0.94,  and  n  =  1/5. 


A  simple  scaling  law  for  the  post  collapse  behavior  of  the  wake  is  thus 
obtained  from  the  two  analyses  when  the  vertical  size  of  the  wake  is  divided 
by  the  Froude  number,  F  ,  which  is  proportional  to  the  ratio  of  the  turbulence 
energy  to  the  potential  energy.  In  physical  terms,  at  a  given  wake  position 
X/D,  the  vertical  extent  of  the  wake,  ~£v  / D  ,  which  is  directly  proportional 
to  the  Froude  number,  F  ,  increases  as  the  wake  turbulence  increases  with 
submarine  velocity  u.  and  decreases  as  the  potential  energy  (D/T)  increases. 
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s  In  the  use  of  Equation  (9)  just  described,  there  is  some  uncertainty  about  the 

values  for  K  and  rt  while  the  degree  of  mixing  in  the  wake  ft/ CL  needs  to 
be  investigated  further. 

In  Figure  12,  a  comparison  is  made  between  predictions  of  thu 
vertical  dimensions  of  the  wake  using  Ko's^  computer  analysis  and  the  power 
law  growth  rates  developed  on  the  present  program.  The  vertical  size 
is  plotted  against  wake  position  X/D  for  F  =  250.  The  Froude  number  F 
is  inversely  proportional  to  the  square  root  of  the  Richardson  number.  For 
typical  ocean  conditions  and  submarine  dimensions,  a  value  of  250  for  F 
corresponds  to  a  submarine  speed  of  roughly  five  knots. 

The  curve  for  unstratified  flow,  0jD  =  K  (X/D)  plotted  in 
Figure  12  gives  the  extent  of  the  wake  up  to  collapse.  A  family  of  curves 
obtained  from  Figure  11  for  different  degrees  of  mixing  ft  / CL  gives  the  size 
of  the  wake  after  collapse.  Comparing  the  ft / CL  =0.94  curve  with  Ko's^ 
analysis,  shown  by  the  dotted  line  in  Figure  12,  excellent  agreement  is  ob¬ 
tained  except  in  the  immediate  vicinity  of  collapse.  At  the  intersection  of  the 
unstratified  flow  and  post  collapse  curves,  a  maximum  difference  of  roughly 
12%  is  observed.  Clearly,  a  fairing  of  the  two  power-law  curves  would  give 
closer  agreement  with  Ko's  result  in  this  region.  The  important  point  is  that 
both  analyses  give  the  same  decay  rate  of  vertical  size  with  distance  in  the 
wake . 

To  estimate  the  size  of  a  submarine  wake,  it  is  convenient  to 
express  the  dimensions  in  terms  of  the  time  after  submarine  passage. 

Before  collapse,  the  growth  of  the  circular  wake  given  by  «  K  (X/B) 

can  be  expressed  in  the  form 

1  1  (10) 

Combination  of  Equation  (10)  with  the  analytical  estimate  of  the  vertical  extent 
of  the  wake  during  collapse  given  by  Equation  (5)  and  confirmed  by  the  experi¬ 
mental  data  in  Figure  10  yields  an  expression  for  4 v  / D  in  terms  of  t/T  . 
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t  saywga^*^ C93?^??gg^g 


jfv_  _  ^nK(f;H  (A 

D  "  ( i  -  /3/ct;Vs  '  T  ' 

(i  : 

Equation  'll)  gives  the  vertical  wake  height  normali.-.'.d  by  the  initial  wake 
diameter  as  a  function  of  the  initial  Froude  number,  tne  degree  of  mixing  xn 
the  wake,  two  parameters  determined  from  the  wake  growth  rate  in  un¬ 
stratified  flow,  and  the  time  after  generation  normalized  by  the  ambient 
Brunt- Vaisala  period.  Substitution  of  the  experimental  values  of  K  =  1.3, 
n  =  1/4  (Figure  2)  and  /3/a~0.93  (Figure  10)  into  Equation  (11)  yields 


\  -  0.5  (F)%  (f  )'V" 


(12) 


4.5.2  Vertical  Asymptote  and  Maximum  Extent 

The  vertical  size  of  the  wake  can  be  estimated  using  Equation  (12) 
between  collapse  and  t/T  «  1.0  at  which  point  measurements  such  as  those 
in  Figures  8  and  10  indicate  that  the  wake  vertical  size  may  reach  an  asymptote 
vp  .  Putting  t/T  =  1.0  in  Equation  (12)  gives 


*  0.5  ( F) 


v4 


(13) 


Equation  (13),  which  requires  further  verification,  gives  a  relation  for  the 
asymptotic  vertical  wake  size  4p  in  terms  of  the  initial  size  and  Froude 
number.  For  the  conditions  of  the  present  experiment,  a  =  1.1  cm/sec, 

D  =  0.6  cm,  and  T  =  55  sec,  F  =  ~  =  100  .  As  shown  in  Figure  13, 
Equation  (13)  gives  an  asymptotic  wake  height  ^/D  of  1.6  compared  with  a 
measured  value  of  1.7  in  Figure  8.  For  Schooley  and  Stewart's^  self- 
propelled  model  experiments  in  which  U.  =  45  cm/ sec,  D  =  2.2  cm,  and 
T  =  2.3  sec,  Equation  (13)  gives  a  value  of  1-4  for  -ftp  / D  equal  to  that 
measured. 

A  similar  relation  can  be  obtained  for  the  maximum  vertical  extent 
of  the  wake  by  substituting  t C/T  «  1/3  from  Figure  4  into 
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Equation  (10)  and  noting  from  Figure  10  that  $rn  /  $ os  ~  to  give 

~^rm  /  p  \  '4 

— '  ^  0-9  (F) 

D  (14) 

For  the  present  experiment,  Equation  (14)  as  shown  in  Figure  13  gives 

/D  »«  2 .8  compared  with  a  measured  value  of  2.9  while  for  Schooley 
and  Stewart's  tests  the  measured  and  estimated  values  are  both  2.4. 

Ko's^  analysis  is  shown  in  Figure  13  to  give  the  same  variation  of 
ir^  / D  with  Froude  number,  F  ,  although  the  values  of  1 rm/D  are  slightly 
lower  than  those  measured.  In  addition,  Ko's  analysis  predicts  values  of 

/D  less  than  those  measured  in  the  region  where  an  asymptote  appears 
to  be  reached.  This  couid  account  in  part  for  Ko's  analysis  overpredicting  the 
horizontal  extent  of  the  wake. 

Van  de  Watering  et  al^  report  measurements  of  $ m  /  D  and  ^  /B 

in  terms  of  the  initial  wake  growth  rate.  Use  of  Ko's  assumptions  that 

u‘/ u.  =  0.25  initially  for  Van  de  Watering's  measurements  gives  values  of 
(F)  between  2  and  10,  which  would  correspond  to  submarine  velocities  of  less 
than  one  knot.  Van  de  Watering's  measurements  also  yield  the  (F)  ^ 
variation  given  in  Figure  11  in  which  his  measurements  at  the  highest  value 
of  F  =  10  are  shown  to  agree  well  with  Equations  (13)  and  (14). 

4.5.3  Horizontal  Width 

As  the  wake  collapses  vertically,  the  horizontal  growth  rate  is 

enhanced  so  that  the  horizontal  width  ^  becomes  greater  than  the  diameter 

in  unstratified  flow  .  The  ratio  1 ^  from  the  present  as  well  as  other 

6-11  ^ 

investigations  is  correlated  quite  well  in  Figure  10  by  the  parameter  t/T 
and  reaches  a  value  of  about  1.6  at  t/T  =  2  compared  with  the  vertical  size 
$-v  /  £«}  ~  •  In  other  words,  the  wake  width  at  this  point  is  roughly  five 

times  the  wake  height.  Ko's  theory  is  shown  in  Figure  10  to  overpredict  the 
horizontal  wake  size  for  t/T  greater  than  about  0.7.  This  is  probably  due 
to  Ko's  assumption  of  an  elliptical  wake,  although  the  presence  of  the  sidewalls 
ma,  ave  restricted  the  horizontal  growth  in  some  of  the  previous  experiments^. 
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To  predict  the  horizontal  extent  of  a  wake  during  collapse,  use  can  be 
made  of  the  correlation  of  measurements  in  Figure  10-  The  data  can  be 
fitted  by 


This  expression  when  combined  with  the  unstratified  growth  rate,  Equation  (10), 
gives  for  K  =  1.3  and  ti  =1/4 


Equation  (16)  gives  the  horizontal  enhancement  of  the  wake  growth 
rate  during  the  period  when  the  wake  is  collapsing  vertically.  If,  as  postulated, 
a  vertical  asymptote  is  reached  at  about  t  / T  £  1.0,  it  is  reasonable  to  expect 
that  the  horizontal  growth  rate  would  be  reduced.  In  this  region,  the  wake  width 
$c  is  assumed  to  increase  at  the  same  rate  as  in  unstratified  flow  with 
n  =1/4  but  with  K  =  2.0  because  Equation  (16)  must  be  matched  at 
t/T  =  1.0.  This  gives 


ic 

D 


z.o(F)  +(-f) 


(1?) 


5.  ESTIMATES  OF  SUBMARINE  WAKE  DIMENSIONS 


In  the  past,  wake  growth  rates  measured  in  ballistic  ranges  and  in 
small  scale  laboratory  investigations  of  unstratified  flow  over  fixed  models 
have  been  used  successfully  to  predict  the  wake  growth  for  full-scale  unpowered 
reentry  vehicles.  In  a  similar  manner,  the  laboratory  measurements  of  grid 
generated  and  model  wakes  can  be  used  to  predict  the  wake  growth  behind  a 
self-propelled  submarine  vehicle  using  Equation  (10)  up  .o  collapse 
(  tc/T  ^  1/3  )  and  Equations  (12)  -  (17)  downstream  of  collapse.  With  the 
best  values  presently  available  for  the  uncertain  parameters  used  in  those 
equations  (K  =  1.3,  n  =  1/4,  ft>/(X  =  0.93),  it  is  possible  to  estimate  sub¬ 
marine  wake  dimensions  for  different  speeds,  initial  vake  sizes,  and  ambient 
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stratification.  In  addition,  the  effect  of  varying  the  parameters  K,  rt  ,  and 
/3/<X  determined  from  theory  and  experiment  could  be  calculated. 

Referring  to  Equations  ( 12)  and  (16),  .  it  is  observed  that  bbth  the 

vertical,  »  and  horizontal,  ,  dimensions  of  the  wake  vary  weakly 

*/♦  ' 

with  the  submarine  velocity  a  being  proportional  to  a  .  As  a  conse¬ 
quence,  measurements  of  wake  dimensions  at  equal  times  after  submarine 
passage  would  riot  change  very  much  with  submarine  velocity.  On  the  other  hand, 
the  vertical  extent  of  the  wake  during  collapse  is  very  sensitive  to  the  ambient 

stratification  being  directly  proportional  to  thb  Brunt-Vai&ala  period  'T  , 

1/4 

although  the  asymptotic  height  only  varies  with  T 

The  scaling  relations  established  in  Section  4  from  experimental  data 
for  wake  growth  and  collaps^  reveal  that  the  important  pararrieters  are  the 
Froude  number,  F  ,  and  the  ratio  t  /T  of  the  residence  time  in  the  wake  to 

I 

the  Brunt-Vaisala  period.  This  is  demonstrated  in  Figure  14  where  the  1 

wake  dimensions  #  /D  divided  by  (F),/4  are  plotted  against  '  t  /  T  ,  .  In  this 

manner,  two  unique  curves  are  obtained  for  estimating  the  horizontal  width  and 

vertical  height  of  a  wake  in  stratified  flow.  Verification  of  these  curves  or1 

similar  ones  possibly  modified  by  different  values  for  the  constants  arid  ex-  1 

ponents  must  rest  ultimately  on  comparison  with  submarine  data. 

1  ' 

From  the  equations  and  curves  for  wake'  growth  in  Figure  14,  unique 
curves  can  also  be  determined  for  the  wake  area.  As  shown  in  Figure  15, 
this  is  accomplished  by  plotting  (?)  *  »  the  wake  area  normalized  by 

the  initial  area  and  divided  by  the  square  root  of  the  Froude  number,  against 
the  time  after  wake  generation.  Tpe  curves  show  that  when  the  flow  is  stratified, 
the  inhibition  of  vertical  entrainment  after  collapses  reduces  the  wake  area  to 

approximately  one  half  the  area  anticipated  for  unstratiiied  flow.  The  assump- 

1  1 

tion  of  an  elliptical  wake  during  collapse  is  seen  in  Figure  15  to  lead  to  an  un¬ 
realistic  decrease  in  wake  area  with  time.  This  can  be  avoided  by  assuming 
that  the  wake  becomes  nearly  rectangular  in  shape  soon  after  collapse. 

*  (  *  I 

As  shown  in  Figure  16,  the  equations  developed  (Figure  14)  can  be 
used  for  preliminary  estirhates  of  the  growth  of  the  wake  behind  a  submarine 
travelling  through  a  stratified  ocean.  Parameters  chosen  for  the  calculation 
are  a  submarine  velocity  U  of  12  knots',  an  initial  \Vake  size  D  of  25  feet 
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and  a  Brunt- Vaisala  period  T  of  12  minutes.  Referring  to  Figure  16,  the 
wake  grows  initially  in  a  conical  manner  and  reached  a  diameter  of  about 
100  feet  at  two  minutes.  A  maximum  vertical  height  of  about  110  feet  is 
reached  between  three  and  five  minutes  followed  by  a  collapse  to  the  asymp¬ 
totic  vertical  height  of  approximately  60  feet  at  fifteen  minutes.  During  this 
collapse  phase,  the  horizontal  growth  rate  is  enhanced  so  that  by  fifteen  min¬ 
utes  the  wake  is  roughly,  260  feet  wide.  By  forty  minutes  lifter  wake  generation, 
which  corresponds  to  a  position  nine  miles  downstream  of  the  submarine,  the 
wake  is  estimated  to  be  approximately  300  feet  wide  and  60  feet  high.  From 
Figure  3,  the  turbulent  velocity  in  the  wake  at  this  time  is  estimated  to  be  of 
order  0.5  cm/ sec  so  that  subsequent  dispersion' of  the  wake  constituents 
would  depend  on  the  background  ocean  turbulence. 

1  Estimated,  extremes  in  wake  dimensions  can  be  calculated  from  the 

equations  and,  curves  developed  in  Figure  14.  The  horizontal  and  vertical 
sizes  forty  minutes  after  submarine  passage  are  plotted  in  Figure  17  against 
submarine  velocity  for  different  oceanic  stratification.  The  initial  wake 
diameter  i,s  taken  as  twenty-five  feet.  For  a  submarine  moving  very  slowly 
(  a  =  2  knots)  through  a  highly  stratified  ocean  (T  =  4  min.),  the  wake  height 
would  be  30  feel.  In  contrast,  a  fast  submarine  (  u  =  25  knots  in  a  slightly 

«  i 

sUutified  ocean  'T  ,=  40  min.)  would  produce  a  wake  100  feet  high.  A  wake 
500  feet  wide  would  be  produced  for  u  =  25  knots  and  T  =  4  min.  whereas  the 
width  would  be  only  150  feet  for  U  =  2  knots  and  T  =  40  min. 

With  an  estimate  of  the  size  and  shape  of  the  v  ake  available,  prelim- 

•  I, 

inary  predictions  of  the  concentration  of  a  passive  tracer  introduced  at  the 
submarine  could  be  made.  In  a  given  cross- section,  it  would  be  reasonable 
to  anticipate  maximum  concentrations  at  the  wake  center  with  some  form  of 
exponential  decay  in  the  horizontal  direction  to  zero  at  the  wake  edge  and 
little  variation  in  the  vertical  direction.  Based  on  the  estimates  of  wake  area, 
high  values  of  tracer  concentration  would  be  observed  at  large  distances  from 
the  submarine. 
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e  3  DECAY  OF  TURBULENCE  INTENSITY  (UNSTRATIFIED  FLOW) 
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(a)  AMBIENT  PROFILE 


(b)  t/T  *  0.06  PROBE  FIXED  NEAR 
TOP  WAKE  EDGE 


(e)  WAKE  PROFILE  t/T=*  0.6 


(f)  WAKE  PROFILE  t/T  =*  1.5 


Figure  6  MEASURED  TEMPERATURE  PROFILES  T~  16 
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PRESENT  ANALYSIS  n  =  1/4,-  =  - iTo 

fa*  t/T  (1-/?/*)1/2 


- USING  KO's  THEORY  B  /a.  =  0.9 

•,0  SCHOOLEY,  STEWART  SELF  PROPELLED  T  =  2.8 
SCHOOLEY  2D  QUIESCENT  T  =  7 
SUNDARAM  2D  QUIESCENT  T  =  12 
4  VAN  DE  WATERING  2D  QUIESCENT  T  =  25 
+  PRESENT  EXPERIMENTS  T  =  55 
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TIME  AFTER  WAKE  GENERATION,  t/T 


Figure  10  CORRELATION  OF  WAKE  DIMENSIONS  DURING  COLLAPSE 
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160  240 

WAKE  POSITION,  X/D 


Figure  12  PREDICTIONS  OF  VERTICAL  WAKE  SIZE,  F  =  250 


SCALING  LAWS  FOR  WAKE  GROWTH  AND  COLLAPSE 


GENERATION,  t/T 
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